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ABSTRACT 

A  variable  conductance  heat  pipe  was  designed  and  con- 
structed to  study  the  performance  characteristics  of  the 
gas  loaded  heat  pipe  with  particular  emphasis  on  the  effect 
of  gravity.   Heat  inputs  were  varied  from  10  to  50  watts 
and  the  operating  position  of  the  heat  pipe  was  varied  from 
horizontal  to  vertical  with  the  condenser  above  the  evapora- 
tor.  Condenser  temperature  profiles  are  presented  for 
various  heat  inputs  and  positions  for  the  conventional  and 
gas  loaded  heat  pipe.   The  possibility  of  stratification  in 
a  horizontal  heat  pipe  operating  with  a  non-condensible 
gas  that  has  a  larger  molecular  weight  than  the  working 
fluid  is  discussed  based  on  observations  presented. 

Additionally,  a  one  dimensional  analysis  was  formulated 
and  used  to  compare  with  the  experimental  results.   The 
results  show  a  good  agreement  between  theory  and  experiment, 
indicating  that  the  developed  mathematical  model  adequately 
describes  the  physical  situation. 
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I.   INTRODUCTION 

A.  BACKGROUND 

The  heat  pipe  is  a  device  capable  of  transferring  large 
amounts  of  heat  while  operating  at  a  nearly  constant  temper- 
ature.  The  major  development  of  this  device  has  occurred 
within  the  last  decade  although  the  original  concept  was 
proposed  in  19^.   This  proposal,  by  R.  S.  Gaugler  [Ref.  1] 
of  the  General  Motors  Corporation,  was  not  acted  upon  at 
that  time.   A  similar  proposal  by  L.  Trefethen  [Ref.  2]  was 
made  to  the  General  Electric  Company  in  1962  but  was  not 
investigated.   The  advent  of  the  space  program  led  to  a 
renewed  interest  in  the  heat  pipe  in  1964  by  G.  M.  Grover 
et  al  [Ref.  3]  of  the  Los  Alamos  Scientific  Laboratory. 
The  possibility  of  developing  a  reliable,  efficient  and 
lightweight  heat  transfer  device  for  application  in  space 
equipment  was  too  appealing  to  allow  the  concept  to  again 
be  overlooked. 

B.  HEAT  PIPE  OPERATION 

The  sealed  heat  pipe  is  lined  with  a  capillary  wick 
structure.   The  vessel  is  evacuated  and  partially  filled 
with  a  fluid  referred  to  as  the  working  fluid.   During 
operation  the  working  fluid  is  continuously  evaporating  in 
one  section  of  the  pipe  and  condensing  throughout  another 
section  (see  Figure  1).   The  condensate  is  returned  to 
the  evaporator  section  by  means  of  capillary  force  acting 


through  the  wick  structure.   This  change  of  phase  process 
results  in  a  nearly  constant  temperature  throughout  the 
heat  pipe.   The  near  isothermal  operation  of  the  heat  pipe 
was  reported  in  1965  by  K.F.  Bainton  [Ref.  *i].   The  iso- 
thermal characteristic  of  the  heat  pipe  allows  for  a  more 
efficient  heat  transfer  device  than  that  achieved  with  an 
equivalent  solid  metal  rod.   An  effective  thermal  conduc- 
tivity for  a  lithium  heat  pipe  of  1,000  to  10,000  times  that 
of  an  equivalent  metal  rod  was  reported  by  workers  at  RCA 
in  1967  [Ref.  5L 


wall 


wick 


Heat  pipe  vapor  flow  diagram 
figure  1 


The  introduction  of  a  non-condensable  gas  along  with 
the  working  fluid  alters  the  performance  of  the  heat  pipe 
appreciably.   During  operation  the  non-condensable  gas  is 
driven  to  the  condenser  along  with  the  working  vapor. 
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Unlike  the  vapor  however  the  gas  does  not  condense  and 
thereby  forms  a  high  concentration  of  the  non-condensable 
gas  in  a  portion  of  the  condenser.   As  the  concentration  of 
the  non-condensable  gas  increases  the  amount  of  working 
fluid  present  and  the  condensation  rate  decrease.   The 
temperature  of  the  pipe  in  this  portion  of  the  condenser 
will  therefore  become  less  than  in  a  portion  where  the  non- 
condensable  gas  is  absent  or  is  a  very  small  amount  of  the 
concentration.   Since  the  pressure  of  a  saturated  vapor  is 
directly  proportional  to  its  temperature,  the  amount  of 
working  vapor  present  will  decrease  with  the  temperature. 
For  a  large  enough  concentration  of  non-condensable  gas 
in  a  portion  of  the  condenser,  ambient  temperature  will 
be  approached.   The  small  amount  of  working  fluid  present 
forming  the  partial  pressure  of  methanol  will  not  condense, 
thereby  reducing  the  heat  transfer  area  and  the  conductance 
of  the  heat  pipe. 

If  the  evaporator  temperature  increases,  the  partial 
pressure  of  the  working  fluid  increases  thereby  compressing 
the  non-condensable  gas  into  a  smaller  volume  which  in 
turn  increases  the  heat  transfer  area.   The  pipe's  conduc- 
tance is  thereby  increased.   This  characteristic  of  the 
gas  loaded  heat  pipe  enables  it  to  have  the  potential  to 
experience  smaller  operating  temperature  variations  than 
a  conventional  heat  pipe  for  similar  changes  in  heat  input. 
Therefore,  the  process  by  which  the  non-condensable  gas 
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collects  in  the  condenser  has  a  direct  beaming  on  the 
heat  pipe's  performance. 

Two  analyses  have  been  developed  to  give  insight  into 
the  behavior  of  the  gas  loaded  heat  pipe.   The  first,  the 
flat  front  theory,  was  developed  by  W.  Bienert  [Ref.  6]. 
This  theory  assumes  that  during  steady  state  operation  of 
the  gas  loaded  heat  pipe  the  vapor  and  gas  are  separated 
in  the  condenser  by  a  vertical  flat  front.   The  axial  heat 
conduction  of  the  pipe  is  considered  negligible  in  this 
theory.   The  possibility  of  a  pressure  grac'ient  being 
present  in  the  pipe  is  excluded  and  diffusion  is  not  con- 
sidered.  With  all  of  these  assumptions  an  analysis  of  the 
heat  pipe  performance  is  straight  forward  amd  proves  ade- 
quate in  the  evaluation  of  several  heat  pip°e  parameters. 
Specifically  the  volume  of  the  gas  reservoir  and  the  ratio 
of  condenser  to  reservoir  volume  can  be  determined.   Figure 
2  represents  the  flat  front  theory. 


Vapor-gas  flat  front  model 

figure  2 
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More  recently,  Marcus  and  Edwards  [Ref.  7]  have  developed 
the  diffuse  front  theory.   In  this  analysis  the  contribution 
of  binary  diffusion  between  the  non-condensable  gas  and  the 
working  vapor  Is  considered  along  with  axial  conduction 
in  the  pipe  wall.   In  addition  to  heat  transfer  from  the 
condenser  by  convection  the  contribution  due  to  radiation 
is  included.   The  inclusion  of  these  phenomenon  results 
in  the  vapor-gas  interface  spreading  out  into  a  more 
realistic  profile.   This  analysis  develops  accurate 
temperature  and  concentration  profiles  for  the  condenser 
section. 

This  analysis  was  used  to  formulate  the  TRW  computer 
program  which  is  capable  of  calculating,  for  various  gas 
loaded  heat  pipe  geometries,  among  other  parameters: 

(a)  the  wall  temperature  profile; 

(b)  heat  and  mass  transfer; 

(c)  gas  reservoir  size. 

C.   OBJECTIVE 

When  a  container  is  filled  with  two  vapors  of  dissimilar 
molecular  weight  the  vapor  of  larger  molecular  weight  will 
settle  to  the  bottom  portion  of  the  container  due  to  gravity. 
In  the  operating  gas  loaded  heat  pipe  the  non-condensable 
gas  tends  to  collect  at  the  end  of  the  heat  pipe  away  from 
the  evaporator,  having  been  convected  there  along  with  the 
working  vapor.   An  inherently  unstable  situation  will 
therefore  be  present  in  the  gas  loaded  heat  pipe  if  the 
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condenser  is  above  the  evaporator  and  the  molecular  weight 
of  the  non-condensable  gas  exceeds  that  of  the  working 
vapor.   The  tendency  of  the  gravitational  force  is  to 
oppose  the  convection  of  the  non-condensable  gas,  thereby 
spreading  out  the  concentration  of  the  gas.   This  decrease 
in  the  gradient  of  the  working  vapor  concentration  is 
accompanied  by  a  decrease  in  the  temperature  gradient. 
Therefore,  it  was  the  objective  of  this  work  to  Investigate 
the  operation  of  a  gas  loaded  heat  pipe  with  particular 
emphasis  on  the  effects  of  gravity. 

The  heat  pipe  used  was  designed  to  operate  horizontally, 
vertically  with  the  condenser  above  the  evaporator  or  any 
angle  in  between  these  two  positions.   This  variation  of 
the  angle  of  operation  allowed  a  variation  in  the  axial 
component  of  the  gravitational  force.   The  results  presented 
herein  are  the  product  of  both  data  obtained  from  a  heat 
pipe  designed  and  constructed  for  this  purpose  and  a  one 
dimensional  analysis. 
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II.   DESCRIPTION  OF  EXPERIMENT 

A.   OVERALL  DESIGN 

As  mentioned  earlier  the  axial  conduction  of  the  heat 
pipe  wall  tends  to  spread  out  the  gas  vapor  interface.   To 
reduce  this,  the  heat  pipe  was  constructed  from  a  stainless 
steel  pipe  with  a  wall  thickness  of  0.0635  inches.   Fins  in 
the  form  of  Acme  threads  were  cut  on  a  lathe  along  the 
condenser  section  to  a  depth  of  0.0255  inches. 

One  of  the  two  end  plates  was  soldered  in  place  while 
the  other  was  designed  to  be  easily  removed  from  the  condenser 
end  of  the  heat  pipe.   The  removable  end  plate,  as  shown  in 
Figure  3,  was  in  the  form  of  an  internally  threaded  cap.   The  cap 
was  fastened  to  a  threaded  flange  which  was  welded  to  the  end 
of  the  heat  pipe.   An  0  ring  was  placed  in  the  flange  to 
effect  a  pressure  and  vacuum  seal.   The  condenser  end  cap  was 
designed  with  a  0.25  inch  diameter  tube  soldered  to  its 
external  side.   Through  a  T  connection  the  tube  provided  the 
necessary  inlet  for  a  pressure  transducer,  vacuum  valve  and 
means  of  evacuating  and  filling  the  heat  pipe. 

A  mesh  stainless  steel  screen  around  the  inside  wall  of 
the  pipe  was  used  for  the  wick.   To  retain  the  wick  firmly 
against  the  wall  a  longitudinal  stainless  steel  spring  was 
inserted  along  the  entire  length  of  the  pipe. 

The  evaporator  heating  device  was  a  resistance  heater 
made  of  .125  inch  wide,  flat  Nichrome  strip  powered  by  a 
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direct  current  power  supply.   Figure  *J  shows  the  overall 
design  of  the  heat  pipe. 


Outside 
thread 


O  ring 


End  cap  design 
figure  3 
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B.   DESIGN  CONSIDERATIONS 

The  performance  of  a  heat  pipe  is  affected  by  any 
chemical  reaction  or  decomposition  of  the  wor-king  fluid  or 
the  corrosion  or  erosion  of  the  wick,  springy  or  heat  pipe. 
An  example  is  the  reaction  of  water  and  stainless  steel 
resulting  in  the  release  of  hydrogen  gas.   In  designing  a 
heat  pipe  the  materials  used  for  the  pipe  anc?.  the  wick  must 
be  chosen  to  be  compatible  with  each  other  amd  the  working 
fluid.   For  this  reason,  the  screen  mesh,  pipe  and  spring 
wire  were  all  made  of  the  same  material,  stainless  steel. 

Because  the  heat  pipe  operation  is  dependent  on  capillary 
pumping  the  working  fluids  surface  tension  had  to  be  consid- 
ered.  In  addition,  the  gravitational  force  and  flow  losses 
made  the  consideration  of  the  working  fluid'::;  density  and 
viscosity  important.   Methanol,  whose  properties  are  slightly 
less  desirable  than  water  was  chosen  because  of  its  good 
wetting  characteristic  and  high  compatibility  with  stainless 
steel.   Ammonia  was  rejected  because  of  its  toxicity. 

Following  Marcus  [Ref.  7],  the  screen  mesh  size  and 
dimensions  of  the  pipe  for  a  given  working  fluid  can  be 
related  by  a  pressure  force  balance  equatior.  .   This  equation 
is  based  on  the  principle  that  the  net  capillary  head  in  the 
wick  must  be  greater  than  the  losses  due  to  pressure  drop  in 
the  vapor  and  liquid  plus  any  losses  due  to  body  forces. 
This  force  balance  equation  is  expressed  as: 
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net        body        liquid        vapor 
capillary  -   force   >   pressure   -   pressure 
head        head        drop  drop 


APc  "  APb  >   PL  +  APv 


AID        ^° 

where,   AP   =  -r- 


APb  =  2Rw  pL 

L  yT  mT 
APT  = 


^   J  PL  (Rw2  "  Rv2) 


and      mT  = 


L    hfE 


and  the  vapor  side  pressure  drop  is  assumed  to  be  much  less 
than  the  other  pressure  components  and  therefore  not 
considered. 

The  heat  load  was  anticipated  to  be  in  the  range  of  zero 
to  100  watts  with  the  most  probable  loading  between  10  and  50 
watts.   The  stainless  steel  screen  mesh  was  available  in 
various  sizes  up  to  200  wires  per  inch.   Using  this  infor- 
mation and  the  above  force  balance  equation  an  iterative 
design  ensued  for  an  operating  temperature  of  120  to  170  °F. 
Various  screen  wire  spacings,  lengths  of  pipe  and  number  of 
layers  of  mesh  were  combined  with  various  pipe  diameters  in 
seeking  the  proper  design.   The  design  centered  on  providing 
for  a  pipe  from  one  to  two  feet  long  and  with  a  diameter  in 
excess  of  one  inch  if  possible. 
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With  a  margin  of  safety,  the  design  v/as  tentatively  set 
for  a  two-inch  diameter  pipe  of  overall  length  15  inches  with 
5  layers  of  150  mesh  screen.   This  length  allowed  for  a  4 
inch  heater  section  and  a  one  inch  adiabatic  section.   This 
was  the  tentative  design  because  these  parameters  were  for 
a  heat  pipe  not  loaded  with  a  non-condensable  gas.   To 
achieve  an  optimum  design  for  the  gas  loaded  heat  pipe  the 
design  program  developed  by  Marcus  [Ref .  7]  was  employed. 

In  order  to  use  this  program  the  non-condensable  gas  had 
to  be  selected.   The  choice  centered  around  the  inert  gases 
since  they  were  chemically  inactive  and  provided  a  good 
spectrum  of  molecular  weights.   To  assist  in  emphasizing  the 
effect  of  gravity  on  the  vapor-gas  mixture  a  large  difference 
in  the  molecular  weights  of  the  working  fluid  and  the  non- 
condensable  gas  was  needed.   The  choice  of  krypton  allowed 
a  molecular  weight  ratio  of  krypton  to  methanol  of  approxi- 
mately 2.5. 

With  the  heat  pipe  parameters  tentatively  determined, 
Marcus'  computer  program  was  used  to  first  determine  the 
additional  length  of  pipe  needed  for  the  non-condensable  gas 
and  then  to  give  a  performance  prediction.   Two  inches  were 
added  to  the  condenser  to  allow  for  the  proper  performance 
of  the  heat  pipe  with  the  non-condensable  gas.   A  check  of 
the  preceding  force  balance  equation  used  in  the  design 
process  indicated  the  heat  pipe  would  still  operate  satis- 
factorily with  the  extra  2  inches  of  length. 
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C.   INSTRUMENTATION 

The  instrumentation  of  the  heat  pipe  included  both 
temperature  and  pressure  sensors.   Seventeen  thermocouples 
were  placed  in  line  in  the  grooves  of  the  fins  at  an  interval 
of  approximately  0.650  inches  on  the  top  of  the  pipe.   An 
additional  thermocouple  was  placed  inside  the  pipe.   This 
thermocouple,  enclosed  in  a  stainless  steel  sheath,  was 
inserted  2  inches  through  the  evaporator  end  plate.   The 
condenser  thermocouples  were  0.010  inches  in  diameter  copper 
and  const ant an  wires  welded  with  a  bead  at  the  end.   The 
bead  was  made  with  a  Dynatch  thermocouple  welder.   These 
beads  were  attached  to  the  grooves  of  the  condenser  with  a 
Unitek  model  1065  welder.   The  other  end  of  all  the  thermo- 
couple wires  were  led  to  a  Jones  strip  and  from  there  to  a 
Hewlett  Packard  2010C  Data  Acquisition  System. 

Several  layers  of  Johns  Manville  MinK  insulation  were 
wrapped  around  the  heater  to  minimize  heat  losses.   Two 
additional  thermocouples  were  then  placed  across  a  layer  of 
insulation.   By  knowing  the  temperature  difference  across  a 
layer  of  insulation  as  well  as  the  thermal  conductivity  and 
thickness  of  the  insulation  the  heat  loss  from  the  evaporator 
could  be  estimated. 

To  determine  the  power  input  to  the  heater  accurately  a 
known  resistance  was  placed  in  series  with  the  Nichrome 
heater  strip.   By  measuring  the  voltage  across  the  known 
resistor  the  heater  current  was  determined.   This  current 
was  used  with  the  voltage  reading  across  the  heater  to 
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determine  the  power  input  at  the  evaporator  section.   A 
Hewlett  Packard  digital  voltmeter  was  used  for  these 
measurements . 

The  pressure  was  measured  by  a  Pace  KP-15  transducer  and 
recorded  on  a  Hewlett  Packard  7702B  recorder.   Figure  5  shows 
a  schematic  diagram  of  the  heat  pipe  and  its  associated 
equipment.   Included  in  this  diagram  and  in  Figure  6  is  the 
apparatus  used  to  fill  the  heat  pipe  with  the  required 
amount  of  methanol  and  non-condensable  gas.   The  filling  rig 
used  was  a  previously  existing  device  used  by  J.  Woodard 
[Ref.  8]  which  was  modified  to  include  an  inlet  for  the  non- 
condensable  gas  and  an  additional  vacuum  valve. 

D.   CALIBRATION  AND  ERROR  ANALYSIS 

The-  thermocouples  were  calibrated  in  a  constant  tempera- 
ture bath  with  a  platinum  wire  resistance  thermometer  and  a 
Rosemont  commutating  bridge,  while  the  thermocouple  readings 
were  made  on  a  Leeds  and  Northrop  K-3  precision  potentiometer. 
The  thermocouple  readings  were  found  to  vary  from  the  com- 
mutating bridge  temperature  readings  by  no  more  than  ±0.11°F. 
The  data  acquisition  system  was  then  checked  against  the  Leeds 
and  Northrop  potentiometer  and  found  to  have  an  error  bound 
of  ±  0.4°F.   The  resistor  placed  in  series  with  the  evaporator 
heater  was  also  calibrated  on  the  Rosemont  commutating  bridge. 
The  resistance  was  read  as  2.013  ±  0.0002  ohms.   The  Harrison 
pressure  recorder  and  Pace  transducer  were  calibrated  with  a 
Wallace  and  Tiernan  gauge.   The  error  for  the  transducer  and 
recorder  was  ±0.4  psia. 

22 


POWER 

SUPPLY 


VOi 


;;>ly 


VACUUM 
PU 


nf?^ 

r> 

" 

"I. 

1  vl 

It 

X 

Jl 

J] 


_  J 


PRESSURE 
TRANSDUCER 


a 

K 
R 
Y 
P 
T 
0 
N 


DATA  ACQUISITION 
RECORDER 


pressure: 
recorder 


Diagram  of  heat  pipe  and  associated  equipment 

figure  5 
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Photograph  of  heat  pipe  and  associated  equipment 

figure  6 
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An  error  analysis  was  performed  for  all  the  data 
recorded  as  well  as  the  other  parameters  used  in  computations. 
A  summary  of  the  error  bounds  is  shown  in  Appendix  A  along 
with  sample  calculations.   The  data  points  on  the  graphs  were 
drawn  to  include  the  error  bounds  of  the  information  presented 

E.   EXPERIMENTAL  PROCEDURE 

Prior  to  proceeding  with  the  experiment,  the  amount  of 
working  fluid  to  be  used  was  determined.   The  equivalent 
liquid  volume  of  the  working  vapor  for  a  heat  input  of 
50  watts  was  added  to  the  void  volume  of  the  wick  structure. 
The  volume  of  methanol  required  was  approximately  45  milli- 
liters depending  on  the  safety  factor  used.   Operation  of 
the  heat  pipe  with  30,  '10,  and  finally  55  milliliters  of 
methanol  produced  no  discernible  difference  in  the  temperature 
and  pressure  readings.   Since  wick  dry  out,  caused  by  a 
shortage  of  working  fluid  in  the  wick  structure,  can  cause 
a  breakdown  in  the  capillary  pumping  process  the  choice  was 
made  to  use  55  milliliters  of  methanol.   An  excess  of 
working  fluid  would  not  produce  any  detrimental  effects  on 
the  heat  pipe's  performance.   The  krypton  gas  used  in  the 
experiment  was  research  grade  while  the  methanol  working 
fluid  was  spectrophotometric  quality. 

Prior  to  assembling  the  heat  pipe,  the  pipe,  wick,  end 
cap,  and  spring  were  all  scrubbed  with  acetone  and  rinsed 
three  times  with  distilled  water.   The  assembled  heat  pipe 
was  subjected  to  the  same  process  prior  to  attaching  the  end 
plate.   The  procedure  used  in  preparing  the  assembled  heat 
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pipe  for  operation  started  with  evacuation  of  the  heat  pipe 
and  filling  assembly.   The  filling  assembly,  in  addition 
to  supplying  the  working  fluid  and  non-condensable  gas, 
provided  for  the  connection  of  the  heat  pipe  to  the  evacua- 
tion pump.   A  total  of  approximately  30  inches  of  0.25 
inch  inside  diameter  vacuum  hose  connected  the  heat  pipe 
through  the  filling  assembly  to  the  evacuation  pump.   The 
pump  was  able  to  achieve  vacuum,  read  at  the  pump,  of 
2  x  10~  millimeters  of  mercury. 

When  the  evacuation  was  completed  the  heat  pipe  and 
filling  assembly  were  isolated  from  the  evacuation  pump. 
A  measured  amount  of  methanol  was  then  released  into  the 
heat  pipe.   The  pressure  was  recorded  when  the  temperature 
inside  the  pipe  equalled  ambient  temperature  .   The  krypton 
was  next  released  from  its  high  pressure  cylinder  and  allowed 
to  flow  into  the  heat  pipe,  while  monitoring  the  pressure 
and  internal  temperature.   The  amount  of  krypton  present 
could  then  be  calculated  by  using  the  temperature  and  the 
difference  between  the  pressure  reading  and  the  partial 
pressure  of  methanol  with  the  perfect  gas  equation. 
When  the  internal  temperature  of  the  pipe  was  equal  to  the 
ambient  temperature,  it  was  recorded  along  with  the  pressure 
in  the  heat  pipe. 

The  heat  pipe  was  then  isolated  from  the  filling  assembly 
by  means  of  the  vacuum  valve  on  the  stem  of  the  heat  pipe's 
end  cap.   The  pipe  was  then  positioned  for  the  desired  angle 
of  operation  and  the  heat  input  was  applied. 
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The  temperature  and  pressure  readings  were  monitored 
until  a  steady  state  operation  was  detected.   The  steady 
state  condition  was  based  on  the  condenser  temperatures 
rate  of  change  being  less  than  0.5  °F  per  hour. 

The  temperatures  and  pressure  were  recorded  at  steady 
state,  then  the  heat  pipe  was  either  repositioned  or  the 
heat  input  was  altered  to  initiate  a  new  operating  mode. 
This  process  was  continued  until  sufficient  data  had  been 
obtained. 

Shutdown  of  the  pipe  required  only  the  removal  of  the 
heat  input  and  subsequent  release  of  the  internal  fluids 
when  the  heat  pipe  had  achieved  ambient  temperature. 
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III.   EXPERIMENTAL  RESULTS 

A.   CONVENTIONAL  HEAT  PIPE 

The  heat  pipe  was  initially  operated  as  a  conventional 
heat  pipe  with  only  methanol  present.   Heat  inputs  were 
applied  from  10  to  50  watts  in  10  watt  increments  in  each 
of  three  heat  pipe  positions;  horizontal,  condenser  raised 
lJ5°j  and  vertical.   The  data  presented  in  Figures  7,  8,  and 
9  is  representative  of  all  the  condenser  temperature  readings 
for  the  various  heat  inputs  and  positions  mentioned.   As 
these  Figures  show,  the  entire  condenser  operated  in  an 
isothermal  manner  as  expected. 

Figures  10  and  11  show  the  evaporator  temperature  and 
condenser  pressure  for  a  particular  heat  input.   The  data 
points  so  obtained  for  each  of  the  heat  pipe  positions  are 
connected  with  smooth  curves.   The  evaporator  temperatures 
were  used  to  determine  the  saturated  vapor  pressure  of 
methanol  from  the  data  on  methanol  presented  by  J.M.  Smith 
[Ref.  9].   The  pressures  obtained  from  J.  M.  Smith's  table 
were  all  within  ±  5%   of  the  corresponding  pressure  readings 
of  Figure  10.   Therefore,  the  assumption  of  operation  of 
the  heat  pipe  with  a  saturated  working  vapor  is  confirmed. 
The  temperature  readings  used  to  enter  J.  M.  Smith's  table 
were  from  the  evaporator  while  the  pressure  reading  used 
for  comparison  were  taken  from  the  end  of  the  condenser. 
The  fact  that  the  pressures  are  equivalent,  within  the 
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accuracy  of  this  experiment,  and  were  read  at  opposite  ends 
of  the  pipe  indicates  that  there  is  no  pressure  gradient 
throughout  the  conventional  heat  pipe  during  operation. 
Figure  11  indicates  the  heat  input  to  the  heat  pipe 
as  a  function  of  the  evaporator  temperature  for  the  horizon- 
tal and  vertical  positions.   The  variation  of  the  evaporator 
temperature  for  the  changes  in  heat  input  is  shown  to  be 
linear.   The  slope  for  the  vertical  position  in  Figure  11 
is  less  than  that  for  the  horizontal  position.   This  indi- 
cates higher  evaporator  temperatures  for  the  same  heat 
input  in  the  vertical  position  when  compared  to  the  horizon- 
tal position.   An  explanation  for  this  can  be  found  in  the 
heat  transfer  coefficient  of  the  cylinder  which  is  directly 
related  to  the  Nusselt  Number  by: 


.    NNUL  k 
h  =  — = 


where  the  Nusselt  Number  equals: 


horizontal        NNUL  =  .525(Npr  NGr) 


vertical         NNUL  =  .59(Npr  NQr) 


.25 
.25 


where : 


L3  p2  g  6  AT  C 

m  ND   =  E- 

Gr  Pr  . 

u  k 
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where  L  is  the  significant  length  and  equals: 

the  pipe  diameter  for  the  horizontal  position, 
the  overall  pipe  length  for  the  vertical  position. 

Since  the  overall  length  of  the  heat  pipe  is  approximately 
9  times  its  diameter  the  heat  transfer  coefficient  for  the 
horizontal  pipe  exceeds  that  for  the  vertical  case.   For  the 
same  heat  input  the  decreased  heat  transfer  coefficient 
in  the  vertical  position,  when  compared  to  the  horizontal 
position,  requires  a  greater  temperature  difference  between 
the  heat  pipe  and  the  atmosphere.   This  is  reflected  in 
Figure  11. 

B.   GAS  LOADED  HEAT  PIPE 

After  evacuation  the  heat  pipe  was  filled  again  v/ith 
55  milliliters  of  methanol  and  2.5  x  10   poundmass  of 
krypton  gas.   The  heat  pipe  was  operated  over  a  similar 
range  of  heat  inputs  and  the  same  three  positions.   The 
condenser  temperature  profile  for  various  heat  inputs  and 
positions  are  shown  in  Figures  12,  13  and  14 . 

Figures  12  and  13  representing  the  horizontal  and  M50 
position,  show  only  a  nomial  negative  temperature  slope 
from  the  evaporator  end  of  the  condenser.   This  temperature 
profile  indicated  the  concentration  of  working  vapor  was 
decreasing  very  little  in  the  direction  of  decreasing 
temperature.   It  was  anticipated  that  a  larger  temperature 
gradient  would  be  present,  reflecting  a  more  marked 
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segregation  of  the  working  vapor  and  the  gas.   It  was  further 
anticipated  that  this  front  between  the  vapor  and  gas 
would  spread  out  over  more  of  the  condenser  as  the  angle 
of  the  heat  pipe  increased. 

The  vertical  operating  position,  as  shown  in  Figure  14 
had  a  larger  negative  temperature  slope  than  that  found  in 
the  other  two  positions.   The  vertical  position  appears 
to  indicate  the  absence  of  a  sharp  front  but  rather  a 
mixture  of  the  vapor  and  gas.   The  concentration  of  the 
vapor  does  decrease  more  markedly  with  distance  from  the 
evaporator  in  this  position.   This  indicates  that  inspite 
of  the  gravitational  force  and  the  larger  molecular  weight 
of  the  non-condensable  gas,  the  gas  does  tend  to  accumulate 
in  the  condenser  at  the  end  away  from  the  evaporator.   The 
concentration  of  the  non-condensable  gas  at  no  point  excludes 
the  condensation  of  the  working  fluid  in  the  condenser. 
Therefore,  the  gravitational  force  does  effect  the  variable 
conductance  property  of  the  gas  loaded  heat  pipe. 

In  comparing  Figures  7,  8,  and  9  with  Figures  12,  13 
and  14  it  is  apparent  that  the  change  in  temperature  of  the 
condenser  for  a  variation  in  the  heat  load  was  less  for  the 
gas  loaded  heat  pipe.   This  is  also  indicated  by  the  larger 
slope  found  in  Figure  11  when  compared  with  Figure  16.   This 
is  another  feature  of  the  gas  loaded  heat  pipe  due  to  its 
variable  conductance.   As  these  two  figures  indicate,  the 
operating  temperature  of  the  evaporator  is  higher,  for  the 
same  heat  input,  for  the  gas  loaded  heat  pipe.   Figures  10 
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and  15  indicate  that,  as  expected,  the  pressure  is  greater 
in  the  gas  loaded  heat  pipe  compared  to  the  conventional 
heat  pipe  for  the  same  heat  input . 

The  lack  of  a  significant  temperature  variation  in  the 
condenser  for  the  horizontal  position  was  investigated 
further  after  the  data  had  been  collected  for  all  heat 
inputs  and  positions. 

The  heat  pipe,  in  the  horizontal  position  was  rotated 
around  its  longitudinal  axis  in  order  to  allow  the  thermo- 
couples to  take  readings  along  the  bottom  of  the  heat  pipe. 
In  this  position  the  readings  were  uniform  and  on  the 
average  12  °F  lower  than  those  taken  on  the  top  of  the 
heat  pipe.   The  bottom  of  the  pipe  was  expected  to  be  cooler 
than  the  top  due  to  the  natural  convection.   However,  the 
temperature  slope  on  the  top  of  the  pipe  and  the  absence  of 
such  a  slope  on  the  bottom  indicated  that  the  non-condensable 
gas  and  methanol  may  have  stratified.   That  is,  the  krypton 
may  have  settled  on  the  bottom  of  the  entire  length  of  the 
pipe . 

If  the  gas  vapor  had  not  stratified  the  nominal  tempera- 
ture gradient  present  on  the  top  of  the  pipe  would  have 
been  present  on  the  bottom  of  the  pipe.   With  stratifica- 
tion, the  heavier  krypton  would  settle  to  the  bottom  portion 
of  the  horizontal  condenser.   This  would  greatly  reduce 
the  condensation  rate  if  not  completely  stop  it  in  this 
section.   The  bottom  of  the  heat  pipe  being  greater  than 
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ambient  temperature  can  be  accounted  for  by  the  circumferen- 
tial conduction  of  the  pipe. 

To  further  investigate  this,  temperature  sensitive 
liquid  crystals  were  applied  over  a  2  inch  wide  black  paint 
strip  which  had  been  applied  to  the  entire  length  of  one 
side  of  the  condenser.   The  liquid  crystal  has  the  charac- 
teristic of  changing  from  clear  to  a  color  over  a  small 
temperature  range.   The  crystals  used  were  R53.   These 
indicated  three  distinct  temperature  ranges  by  three 
different  colors  as  indicated  below  [Ref.  10]: 

Color  Temperature   °C 
Red  50.5  ±   0.5 

Green  52.1  ±  0.5 

Blue  53.3  ±   0.5 

With  the  liquid  crystals  on  the  pipe  the  heat  input 
was  varied  until  the  general  temperature  range  for  crystal 
coloring  was  achieved.   Figure  17  shows  two  photographs  of 
the  isotherms  observed  for  the  indicated  power  input.   In 
the  case  of  a  heat  pipe  the  isotherm  represents  a  line  of 
constant  working  vapor  density.   The  fact  that  the  isotherm 
lines  run  down  the  side  of  the  pipe  in  a  near  horizontal 
manner  indicates  the  working  vapor  density  varies  in  the 
same  manner.   Stratification  to  some  degree  may  therefore 
be  present. 


A3 


" 


Heat  input  equals  21.70  watts 


Heat  input  equals  23.6  watts 

Photographs  of  liquid  crystals  on  condenser 

figure  17 


W 


The  effect  of  natural  convection  upon  the  slope  of  the 
iostherms  has  not  been  fully  investigated  at  this  time  and 
v/ould  provide  more  insight  into  this  phenomenon.   The 
implications  of  the  possibility  of  stratification  will  be 
discussed  in  Section  V. 
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IV .   ANALYSIS 

A.   THEORY 

The  analysis  developed  in  this  section  is  for*  the  con- 
denser section  of  a  gas  loaded  heat  pipe.   The  development 
is  through  the  application  of  the  continuity  and  momentum 
equations  for  one  dimensional  flow,  the  energy  equation  for 
the  wall,  wick,  and  fin  heat  transfer  as  well  as  Fick's 
law  of  diffusion  and  the  perfect  gas  law.   The  development 
resulted  in  four  first  order  differential  equations  which 
were  solved  simultaneously  with  the  use  of  a  Runge  Kutta 
routine.   The  formulation  of  this  analysis  is  outlined 
below.   It  should  be  noted  that  as  Figure  18  indicates,  the 
origin  of  the  coordinate  system  is  at  the  end  of  the  pipe 
away  from  the  evaporator. 


Heat  Pip© 
figure  18 
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Starting  with  the  integral  form  of  the  momentum  equation 
as  applied  to  the  section  shown  in  Figure  19: 


P,*V 


ev  +  ^sm 

°  ox 


P  +  dP 


figure    19 


(1) 


I    Fy   =    /    pV(V-dA.) 


In  the  following  analysis  the  viscous  term  has  been  left 
out  of  the  above  equation.   This  simplifying  assumption  is 
based  on  the  Reynolds  number  being  very  small  in  the 
operating  range  of  the  heat  pipe.   For  a  typical  heat  load 
of  20  watts: 


V  D 
Rn  =  — -  =  50.0 
D     v 


where,     V  = 


Q 


hfg  Ac  p 


V  = 


68.26 


- 


(481)  jj-  (1/6)"  (1/20) 


V  =  .04  ft/sec 
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For  fully  developed  laminar  flov,7  the  friction  factor  f 
can  be  represented  by  64/R„  .   This  results  In  a  shear  force 

P  =  8uti-V  dx 


Considering  the  entire  one  foot  length  of  condenser 
the  force  is  in  the  range  of  10~  pounds  or  less.   Therefore, 
the  shear  force  can  be  neglected. 

Referring  to  Figure  (19),  equation  (1)  becomes: 


P  -  [P  +  |f  dx]  +  pg  sin  0  dx  =  -pV2  +  [pV  +  |^  (PV)  dx] 

[V  +  ?  dx]  +  m  V   7T  D  dx 
dx  a 


where  P  is  the  pressure  in  the  condenser  section,  V  is  the 

vapor-gas  velocity,  V  is  the  velocity  of  the  vapor,  p  is  the 

a 

vapor-gas  density,  g  is  the  gravitational  acceleration  and 
D  is  the  diameter  of  the  vapor  space. 

Dividing  by  A  dx,  and  taking  the  limit  as  dx  approaches 
zero  gives: 


,~s  dP  ,      ,   a    „  dV  ,  „  d(pV)    •  „  ttD 

(2)   -  -=—  +  pg  sin  6  =  pV  -r—  +  V  — ^ — -  +   m  V  -r— 
dx   K&         K   dx      dx        a  A 

c 


The  continuity  equation  for  the  mixture  is 


pVA   =  m  tt  D  dx  +  [PV  +  ^—  (pV)  dx]  A 


or 


(3)        ra  15  +  !_  (pv)  =  o 

c 
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Applying  Pick's  First  Law  to  species  a 


oo  Pava  =  pav-p-sg 

where  J-J  is  the  diffusion  coefficient  for  the  vapor-gas 

mixture,  po  is  the  density  of  the  vapor  and  w  is  P^/p. 
a  a 

An  alternate  form  of  this  equation  is : 
(5)         V(w-l)  -tfgi 


The  pressure  in  the  heat  pipe  condenser,  P,  is  the  sum 
of  the  partial  pressure  of  the  vapor,  P  ,  and  the  partial 
pressure  of  the  gas,  P.  : 

P  -  Pa  +  Pb 

The  perfect  gas  law  can  be  used  for  the  non-condensable 
gas.   This  is  possible  because  the  operating  temperature  and 
pressure  in  the  heat  pipe  are  well  below  the  critical  point 
of  the  gas.   However,  during  heat  pipe  operation  the  working 
fluid  vapor  is  in  a  saturated  state  which  precludes  use  of 
the  perfect  gas  equation.   Instead  the  pressure  of  the 
working  fluid  vapor  can  be  represented  as  a  function  of 
temperature  only.   Therefore, 

(6)  P  =  f(T)  +  pb  A  T 


(7)  P  =  f(T)  +  =£-  p  (l-w)T 

Kb 


^9 


where  R  is  the  universal  gas  constant,  M,  is  the  molecular 
weight  of  the  gas  and  f(T)  is  an  equation  based  on  the 
saturated  vapor  data  presented  by  J.M.  Smith  [Ref.  9]. 

Consider  the  section  of  wick-wall  shown  in  Figure  20  for 
an  energy  balance. 


0" 


q+dq 


V  //////////  //J  ///  ////  UTT7 


D 


h(T-TcC) 


-wick 


<=  q 


-wall 


figure  20 


Assuming  the  temperature  difference  between  the  wick  on 
the  inside  of  the  pipe  and  the  outside  of  the  pipe  is  negli- 
gible in  the  radial  direction: 

(8)   q  +  m  hf  tt  D  dx  =  q  +  dq  +  hCT-T^)  tt  D   dx 

where,  q  is  the  heat  being  conducted  axially  by  the  pipe, 
m  is  the  rate  of  condensation  of  the  working  fluid  per  unit 
area,  h~   is  the  latent  heat  of  condensation  of  the  working  . 
vapor,  h  is  the  heat  transfer  coefficient,  T^  is  the  ambient 
temperature  and  D  is  the  outside  diameter  of  the  heat  pipe. 
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From  Fourier's   conduction   lav/: 

(9)  q   =   -k  A     ~ 
w/  ^  p   dx 

where  A   is  the  cross  sectional  area  of  the  wick-wall  and 
P 

k  is  the  thermal  conductivity  in'  the  axial  direction  of  the 
wick-wall  combination.   The  evaluation  of  the  thermal 
conductivity,  k,  and  the  heat  transfer  coefficient  are 
developed  at  a  later  point  in  this  section. 

Substituting  equation  (9)  into  equation  (8): 

ri2T   h  As    ™  hfP-  *  D 

(10)  2—i  =  -  (T-T  ^  — - ^ 

^lu;  ,2    k  KX        «'  L  A      k  A^ 

dx  p        p 

The  equations  necessary  for  the  development  of  the 
analysis  are  summarized  below: 


(2)     -  j-  +  pg  sin  8  s  pV  j-  +  V  -~ — -  +  -. —  —  m 
dx    Kto         K   dx      dx     A_  w 

c 


(3)      ™|e  =  _^(pV) 

c 


(5)       v  =  T7-l^-3!ic 


(w  -  1)  °^  d 


x 


(6)    P  =  f(T)  +  g£  (l-w)T 

nb 

d2T      •  hfs    *  D    h  As 

(10)    SL-i  =  _  m  -IS +  -  (T-T  )  — — 

K±U)         ,2      m   k  A      k  ^     <*>;  L  A 

dx            p  I 


and, 

(11)  P  =  - 

K   w 
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where  S  is  an  equation  for  p   as  a  function  of  T,  based  on 
the  data  presented  by  J.M.  Smith  [Ref.  9]. 

As  noted  in  the  above  equations,  the  diffusion  coeffi- 
cient,!^ ,  has  been  taken  as  constant.   The  diffusion 
coefficient  equation  for  low  pressure  is  based  on  the 
following  equation  [Ref.  11]: 

.     1.858  10"d  T1^  [  \        d    ] 

P  c122  fiD 

This  indicates  that £Q  does  vary  somewhat  with  temperature 
However,  the  above  assumption  has  been  made  in  order  to 
simplify  the  analysis.   A  mean  diffusion  coefficient  for  the 
temperature  range  involved  was  used  in  the  program. 

The  following  non-dimensional  variables  were  used  in  the 
above  set  of  equations: 


p* 

= 

P 

P 
0 

• 

P 

= 

P 
po 

* 

P 

= 

w 

* 

V 

= 

V 
Vr 

m 

where , 

Vr 

r 
Po" 

and, 

*r 

= 

Q 

fg      c 
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*    m 

m    =  pr 

r 


*    x 

X   =  L 


f*  = 


Po 


Substituting  these  into  the  above  equations: 


(12)   p   sin  0  -  [j-gj-]  -t  =  [(___)   g-] 

^o&   dx      ttD  p  h„    & 

o  fg 

[PV  <  +  v*  "tpV)] 

dx         dx 

+  [(^-)2^]!^ 

ttD  p  h„     & 
o  fg 


(13)   m  s  -  ir  — ?r  (p  v  ) 

dx 

Cw"1)       4  L  Q     dV 

(15)   P   =  f  (T)  +  [  M  g   °]  P   T  (1-w) 

b  o 

,2*    hLA     *    T       iit2^-* 

(16)  ^2  =  TTT^  (T  "  t^  "       Q  i 

dx  ^     K  p         o    D  k  A   T 
^  p   o 


Substituting  for  the  non-dimensional  coefficients  in  the 

preceding  equations:        p 

o 


%  =  p0  g  L 
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N2  "  cViT-]2  3: 

7:0  Pohfg 


n    =  Ok 

W3  D 


m  hL      s 

N4  =  IT  A" 
P 


N      =      *    Tj2  Q 


5  Dk   Ap  TQ 

7rI)3pohfg<£) 

6  16QL2 

n   P  T 


it  Dp  h_ 


3W  6         m: 


the  equations  become: 


fxv\  *    -       a        xr   dP*    KT   r  V  dV*     *  d(o*V*)1 

(17)  P  sin  0  -  N1  — w   =  N2  [p  V  — r  +  V   — *-» — 3— ^-J 

dx  dx         dx 

(18)  m  -  -. -A.  S-,.  <p  v  ) 

N3  dx 

(20)  P*  =  f*(T)  +  N7(l  +  w)P  T* 

d^T  *     *        •* 

(21)  ^i   -  N,,  <T  -T.  )  -  N  m 

dx  ^  ° 


5^ 


Investigating  the  relative  magnitude  of  the  coefficient 
of  the  momentum  term  in  equation  (17)  with  respect  to  the 
other  coefficients  for  representative  values  of  the  operating 
heat  pipe,  one  finds: 

_  r     4x10x3.2113 n2r     1      -. 

N2    L  470x3.  li|(i7B)(l/31.1)J  L(32.2)(3600)J 

N2  =  2.4xio"5 

A  similar  investigation  of  the  other  terms  indicated 

-5 
that  Np  was  at  least  a  factor  of  10.   times  smaller  than  the 

value  of  the  other  non-dimensional  coefficients.   This  led  to 

the  simplifying  assumption  that  the  momentum  term  could  be 

excluded  from  the  analysis  v;ithout  loss  of  accuracy.   The 

computation  of  this  coefficient  was  retained  in  the  final 

program  in  order  to  serve  as  a  check  on  its  relative  value. 

The  momentum  equation  then  reduces  to  a  force  balance: 


*  dP 

(22)  p  sin  6  -  N-.  ^  =  0 

1  dx 


Differentiating  equation  (20): 


(23)  <  .  dfVi  dr^  +      (i^)(s*  +  <  T*)  < 

dx     dT    dx      '  dT      dx 


1   *  *  dw 

w       dx 
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Subtracting  this  from  equation  (22): 


1   S*     ft    df*  dT*  .  TT   ,1-Ww  *  ,  dS*   ^  dT* 
(2*0   —  —  sin  6  =  — r  — r   +  N7  (——MS   +  — w   1* )         ?r 

Nl  w  dT   dx      '    W         dT       dx 


1   *  *  dv; 
-  N7  ~   S  T   ^V 

'  vT       dx 


Solving  equation  (2*1)  for  — w   and  differentiating 

dx 

dx       7  ax 

2  * 

+  [W  G  +  w(i-w)al.](SII.)2 

N?   3  t>   dx 

r2w   p       ,     n       PvTNP    -,    dT      dw 
+    Lvj—  G-,    4     (i  -  ^v:;Oh  J    — r        ?r 
N7      -1  dx     dx 

+       1       [_  ^  sin   G  ^V  +   -^  sin   e   fflj.] 
K1N7  t  dx  T  dx 


where , 


P  df  1 

x         dT      S    T 

dV  SV  _  <  <  T*  .  <  S* 

dT                         dT      dT  dT 

G^    =    ■ TS 

3  (ST) 


r      -     1    +dS      L 
Gu  -  —*  +  — w  nr 

H        T  dT      S 

d^S  *         S      dS  ,dS    ,2 

To  s      +   IF  5F  -  ( W> 

dT  T      dT  dT 

G5    =  -^ 
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Substituting  equation  (18)  into  equation  (21) 


(26)       ^4_  =  N  (T  _T<w  )  +  ^  ^(p  v  ) 
dx  .         3  dx 


Differentiating  equation  (19)  and  substituting  into 
equation  (26) : 


,07x       d^T   _  M  ,*  „  **    M  M  rdS     1     dw   dT 

(27)      -tf2  =  VT  "T~  }  +  N5N6[~^"  w(w-l)  —*  —* 
dx      4  p  o  dT  wvw   ±;  d   dx 


2  (v:  -  1 )  *,dw  v2  S*  d2w  , 
~2"  —2  S  (— r)  +  w(w_1)  — F2  J 
w  (w  -  1)     dx  dx 


Substituting  equation  (25)  into  equation  (27) 


(28)     o^^-Vw-lXT'-O  +  Ou^Sw, 

dx  dx  dx 


-  G^)2  +  G13(<)2 
"^  dx       ±:i  dx 


„       n  dw   .  „       0  dT 
-  G,|.  sin  0  — j-  +  G..  r  sin  6  — w 

1H  dx     ±p       dx' 


where, 


G1Q  =  [(w  -  1)  -  -g-2-  wS  G-l  +  (w  -  DN5N6S  G^] 

G   =  [_5_1  ^   _56  s*   +      (l-2w)  s*   -, 
11      w   ,T*     N7      1     5  6    w       4 

0     rT  >T    2w  -  1    * 

G12  =  N5N6  —, 77  S 

w  (w  -  1) 
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N,-N 


G 


13 


=  [-£-—  wS  GQ  +  Nrl\V  (w-1)  S*  Gc] 


N 


3    5  6 


N5N6  1  S* 
'1U    N-^  W  ^ 


G 


N5N6  S* 


15 


N1N?  ^2" 


Simplifying  equation  (25): 


(29)     d2w    P    d2T*       r<*T*x2.  r    dT*  dw 

7-W2  "  G20  7^2  +  G21(T~^)  +  G22  "T^  7^ 
dx         dx         dx         ax   dx 


n  n  dw     ~     .   0  dT 

-  G?~  sin  G  — 3-  +  Gph  sin  0   — y 


dx 


dx 


where, 


G20  =  lW~   Gl  +  w(1  "v;: 


'k 


G21  =  [f^G3'  +  ^l-^°53 


G22   =  [IF  Gl  +  (1  -2w)(V 


r   -  r  1       1  1 

b23    LN1N?  ^J 


G24  "  [N^"  ^2] 
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Equations  (28)  and  (29),  along  with  their  associated 
non-dimensional  coefficients,  form  a  set  of  two,  second 
order,  simultaneous  non-linear  differentia]  equations.   As 
mentioned  earlier  these  equations  were  solved  using  the 
IBM  360/67  computer  with  a  Runge  Kutta  routine.   The  entire 
program  was  put  in  double  precision  to  avoid  as  much  machine 
round  off  error  as  possible. 

The  value  of  the  heat  transfer  coefficient  for  natural 

convection,  h,  was  calculated  based  on  the  Musselt  Number, 

N    [Ref.  12]. 
nux 


25 
where,        N    =  m(n  „  N  )' 
9  nux      gr  pr 


where,        m  =  .59  for  the  vertical  position 
and,  m  =  .525  for  the  horizontal  position 

The  Nusselt  Number  for  the  heat  pipe  in  the  horizontal 
position  was  found  to  be  11.59,  while  for  the  vertical 
position  the  Nusselt  Number  equalled  67.84.   These  values 
were  used  in  the  following  equation  to  find  the  natural 
convection  portion  of  the  heat  transfer  coefficient,  h  : 

N    k 
,      nux 

To  this  was  added  the  value  of  the  radiation  component 
of  the  heat  transfer  coefficient,  h  ,  resulting  in  the 
following  final  equation: 
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a  e  (T^   -  T  S 

N  00   ' 

h  =  hc  +    (T-T  ) 


Two  values  for  a  condenser  temperature  of  600° R  and  an 
ambient  temperature  of  5^0° R  are  indicative  of  the  relative 
values  used: 

horizontal  position  h  =  1.4l 
vertical  position    h  =  1.00 

In  determining  the  axial  thermal  conductivity  of  the 
wall  and  saturated  wick  combination,  k,  the  conductivity  of 
the  saturated  wick  was  first  determined  and  then  combined 
with  the  thermal  conductivity  of  the  stainless  steel  wall 
in  a  parallel  manner.   The  determination  of  the  saturated 
wick  thermal  conductivity  could  be  viewed  as  either  a 
parallel  or  series  combination  of  the  wire  screen  and  the 
working  fluid  [Ref.  7].   The  two  values  obtained  are: 

k«*  KL 

=  .296 


series   $k  +  (1-4))  kT 

S  Li 


k    ,,  _  =  (1  -  <J))k  +  <|>kT   =  5.80 
parallel       Y   s    y  L    •* 


where  <{>  is  the  void  fraction  in  the  wick,  and  is  equal 
to  .387. 

In  assuming  the  value  is  that  of  the  parallel  case  and 
combining  this  with  the  thermal  conductivity  of  the  stainless 
steel  wall,  the  largest  possible  value  of  thermal  conductivity 
for  the  wick-wall  combination  is  obtained.   This  is  what  the 
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designing  of  the  fins  was  intended  to  reduce.   The  axial 
conductivity  used  was  therefore: 


k  =  15.2     BTU 


hr  ft  °F 


B.   RESULTS 

The  computer  program  developed  to  solve  the  set  of 
simultaneous,  first  order,  non-linear  differential  equations 
is  included  after  Appendix  A.   Several  solutions  obtained 
from  this  program  are  included  in  Figures  21  and  22.   The 
solutions  v;ere  obtained  for  a  specific  heat  input  by  assuming 
a  temperature  for  the  condenser  at  x  equals  zero.   As  a 
first  estimate  the  condenser  temperature  from  the  experimen- 
tal phase  for  the  same  heat  input  was  used.   The  computer 
solution  indicates  in  the  output  the  amount  of  heat  the  pipe 
can  transfer  at  this  temperature  divided  by  the  amount  of 
heat  specified  as  being  the  input.   The  condenser  temperature 
was  varied  until  this  quotient  equalled  one. 

The  analytic  solutions  shown  in  Figures  21  and  22 
indicate  a  nearly  isothermal  condenser.   The  condenser- 
evaporator  interface  is  only  from  0.01  to  0.2$  greater  in 
temperature  than  the  other  end  of  the  condenser  in  these 
figures. 
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V .   DISCUSSION 

A.   COMPARISON  OF  ANALYSIS  AND  EXPERIMENT 

Figures  12  and  lh   present  condenser  temperature  profiles 
based  on  the  operation  of  the  gas  loaded  heat  pipe.   Figures 
21  and  22  present  the  condenser  temperature  profiles  based 
on  the  analysis  developed  in  the  preceding  section.   The 
positions  and  heat  inputs  are  the  same  for  Figures  12  and  21 
and  the  same  for  Figures  14  and  22. 

The  analytic  solutions  all  have  higher  condenser 
temperatures  than  the  data  obtained  from  the  operating  heat 
pipe  for  identical  heat  input  and  position.   The  analytic 
condenser  temperatures  exceed  the  operational  temperature 
in  °F  by  7.2  to  19  percent.   The  larger  difference  between 
the  analytic  and  operational  condenser  temperature  is 
found  in  the  vertical  operating  position. 

The  heat  transfer  coefficient,  h,  if  increased  would 
reduce  the  condenser  temperature  since  for  a  given  heat 
transfer,  Q: 

Q  =  h(T  -  T  ) 

"■  »  CO 

The  values  of  h  determined  earlier  were  based  on  natural 
convection  over  a  symmetric  cylinder.   The  addition  of  fins 
and  the  removable  end  plate  with  its  associated  flange  on 
the  pipe  served  to  disturb  the  flow  over  the  pipe.   This 
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disturbance  of  the  natural  convection  over  the  heat  pipe 
would  have  more  effect  in  the  vertical  position  where  the 
convection  flow  travels  the  length  of  the  pipe.   Any  such 
disturbance  would  tend  to  allow  greater  mixing  in  the  con- 
vective  flow.   This  would  allow  a  greater  amount  of  heat 
to  be  transferred  at  the  same  condenser  temperature  or  the 
original  amount  of  heat  transferred  at  a  slightly  reduced 
temperature.   This  is  the  effect  of  increasing  the  value 
of  h  in  the  preceding  equation.   Increasing  the  value  of 
the  heat  transfer  coefficient  used  in  the  analytic  solution 
will  reduce  the  condenser  temperature  obtained.   This 
appears  to  therefore  be  one  of  the  reasons  for  the  differ- 
ence in  the  analytic  and  operational  condenser  temperatures. 

The  analytic  condenser  temperature  was  nearly  isothermal 
as  noted  earlier.   This  again  differs  from  the  operational 
condenser  temperatures  obtained.   The  condenser  temperature 
decreases  with  increasing  distance  from  the  evaporator  in 
both  the  analytic  and  operational  solutions.   However,  the 
analytic  solution  did  not  have  as  large  a  decrease  as  the 
operating  heat  pipe  did  for  the  vertical  position. 

The  difference  in  the  temperature  slopes  of  the  analytic 
solution  and  operating  heat  pipe  are  not  of  such  magnitude 
to  require  the  analytic  solution  to  be  discarded.   In  fact, 
both  solutions  indicate  the  absence  of  any  front.   Therefore, 
the  assumptions  of  negligible  viscous  and  momentum  forces 
in  the  analysis  are  justified. 
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The  analysis  is  unable  to  give  any  insight  into  the 
possibility  of  stratification  in  a  horizontal  gas  loaded 
heat  pipe,  except  to  confirm  the  lack  of  any  vapor  gas 
front.   The  possibility  of  stratification  in  an  operating 
heat  pipe  cannot  be  overlooked.   The  design  of  the  gas 
reservoir  may  require  alteration  to  allow  the  stratified 
gas  to  flow  in  and  out  with  heat  input  and  pressure 
changes. 

The  heat  transfer  characteristic  of  a  gas  loaded  heat 
pipe  with  stratification  will  differ  from  the  flat  front 
and  diffuse  models  discussed  earlier.   The  stratified  heat 
pipe  may  produce  some  of  the  characteristics  of  the 
diffuse  front  model  with  the  front  turned  90°  in  the 
pipe. 

B.   RECOMMENDATIONS 

Further  investigation  with  the  computer  program  has 
indicated  a  heat  pipe  of  different  dimensions  from  that 
used  in  this  experiment  may  yield  more  interesting  results. 
A  longer  pipe  of  smaller  diameter  with  a  higher  heat  trans- 
fer coefficient  produces  a  marked  vapor-gas  front  in  the 
analytic  solution.   With  such  a  diffuse  front  the  rotation 
of  the  pipe  through  90°  could  be  attempted  to  observe  the 
gravitational  effect  on  the  vapor-gas  front. 

In  addition,  other  inert  gases  could  be  introduced  into 
the  existing  heat  pipe  in  place  of  krypton.   Lighter  gases 
may  allow  the  vapor-gas  front  to  become  more  apparent  while 
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also  providing  a  comparison  for  the  data  in  this  work 
concerning  the  possibility  of  stratification.   The  use  of 
liquid  crystals  could  continue  to  provide  insight  into 
this  phenomenon. 
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APPENDIX  A 


SAMPLE  CALCULATIONS  OF  ERROR  ANALYSIS 


2 
2     V 
Q  -  I2R  -  \ 


WQ=  [(f  Wv)2+  ("^WR)231/2 


where  V  is  the  voltage,  R  is  the  resistance,  WR  is  the 
error  bound  in  the  resistance  and  W„  is  the  error  bound 
in  the  voltage. 

Dividing  both  sides  by  Q: 


W  W    1/2 


WQ    r,    2      nr,2    ,  .003, 2, 1/2 

T  =  [(icrnra-5  x  -05)    +  (276rf>  ] 


W 

-£  =   0.0093 


W0 

-f  =    0.93   J6 
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Quantity 

T  Temperature 

P  Pressure 

Q  Heat  input 

L  Length 

D  Diameter 

gQ  Diffusion  coefficient 

h  Heat  transfer  coefficient 

g  Gravity 

hf  Latent  heat  of  vaporization 

o  Surface  tension 

p  Density 

v  Kinematic  viscosity 

V  Voltage 

R  Resistance 

0  Void  fraction  of  wick 

m  Condensation  rate 


Error  bound 

±    0.5   F 

±    0.4   psia 

±   0.93$ 

±    0.004    feet 

+    0.0004    Teet 

±  i.Qh% 

±   2.8% 

±    0.01    ft/sec2 

±   0.0  3   BTU/lbm 

±    2    x   10~5    lb/ft 

±0.01   lbm/ft3 

±    3   x   10-i4    ft2/hr 

±    0.05    volts 

±    0.00  3   ohms 

±    1.0? 

±    2.4? 
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A  variable  conductance  heat  pipe  was  designed  and  constructed  to 
Ldy  the  performance  characteristics  of  the  gas  loaded  heat  pipe 
Eh  particular  emphasis  on  the  effect  of  gravity.   Heat  inputs  were 
tied  from  10  to  50  watts  and  the  operating  position  of  the  pipe 
r^s  varied  from  horizontal  to  vertical  with  the  condenser  above  the 
r&porator.   Condenser  temperature  profiles  are  presented  for 
'■•rious  heat  inputs  and  positions  for  the  conventional  and  gas 
..aded  heat  pipe.   The  possibility  of  stratification  in  a 
irizontal  heat  pipe  operating  with  a  non-condensable  gas  that  has 
ilarger  molecular  weight  than  the  working  fluid  is  discussed  based 
)  observations  presented. 

Additionally,  a  one  dimensional  analysis  was  formulated  and  used 
;  compare  with  the  experimental  results.   The  results  show  a  good 
ireement  between  theory  and  experiment,  indicating  that  the 
Iveloped  mathematical  model  adequately  describes  the  physical 
s tuation. 
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